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ABSTRACT 



Accretion shocks have been recognized as important X-ray emission mechanism for pre-main sequence stars. Yet the X-ray properties 
of FUor outbursts, events that are caused by violent accretion, have been given little attention. We have observed the FUor object 
ZCMa during optical outburst and quiescence with Chandra. No significant changes in X-ray brightness and spectral shape are 
found, suggesting that the X-ray emission is of coronal nature. Due to the binary nature of Z CMa the origin of the X-ray source is 
ambiguous. However, the moderate hydrogen column density derived from our data makes it unlikely that the embedded primary star 
is the X-ray source. The secondary star, which is the FUor object, is thus responsible for both the X-ray emission and the presently 
ongoing accretion outburst, which seem however to be unrelated phenomena. The secondary is also known to drive a large outflow 
and jet, that we detect here for the first time in X-rays. The distance of the X-ray emitting outflow source to the central star is higher 
than in jets of low-mass stars. 
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1. Introduction 

Variability is a dominant observational signature in pre-main se- 
quence (pre-MS) stellar evolution. Besides the ubiquitous short- 
duration flares related to magnetic reconnection events, various 
types of long-term outbursts are reported: FUor events are char- 
acterized by optical intensity changes of up to 4 mag and a fad- 
ing phase of decades, EXor events are less extreme and shorter 
(months to few years). Both types of outburst are associated with 
a sudden increase of the accretion rate, such that the disk out- 
shines the central star leading to cha racteristic spectral signa- 
tures (e.g. lHartmann & Kenyon|[l996h . The spectra of FUor ob- 
jects resemble that of F-G supergiants, while EXor objects are 
of later spectral type. The FUor and EXor phenomena are prob- 
ably recurrent, about once every 10"' yrs for the FUors and every 
few years in the case of EXors. Only a minor fraction of pre-MS 
stars is classified as FUor or EXor ([Abraham et al.ll20()4tlHerbig| 
l2008 f). Different mechanisms have been proposed as triggers for 
the outbursts: dynamic al interact ion with a close binary com- 
panion (Bon nell&Bast ien 1992; lReipurth & Asprnll2004 . flier- 
mal insta bility in a disk wit h high accretion from a surrounding 
envelope (iBell & Lin 1994j). and changes in the magnetic field 
configuration (Ivan den Ancker et al.ll2004l) . 

In recent years accretion has been recognized to make an 
important contribution to the X-ray emiss ion of pre-MS stars 
dKastner et al.ll2002l:lsterzer & Schmittl200 4). making FUor and 
EXor objects interesting targets for X-ray studies. Plasma tem- 
peratures of up to a few 10^ K can be produced in the accre- 
tion shocks that form when matter is funneled along the mag- 
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netic field lines onto the ste llar surface. In the fir st dedicated X- 
ray survey of FUor objects [Skinner et"an (l2009l) have detected 
two of four stars with XMM-Newton. None of their targets were 
in a state of recent optical outburst during the X-ray observa- 
tion. Only two EXors have been observed in X-rays during an 
optical outburst, with contradictory results dKastner et alj|2006c 
Audard et al. 2005). 

The X-ray spectrum of the prototype FU Ori shows surpris- 
ingly a complex absorption pattern: While the harder emission 
is strongly absorbed and can be ascribed to an embedded stellar 
corona, the origin of the weakly absorbed soft emission is un- 
clear (Ski nner et al . 2006). Possible scenarios include the over- 
laid effect o f a binary companion , accretion shocks, and shocked 
jets. Indeed. lGiinther et alj (l2009l) showed recently that the unre- 
solved soft component in the X-ray spectrum of the pre-MS star 
DG Tau can be explained as emission from a post-shock cool- 
ing zone in the innermost part of its optical outflow. Soft emis- 
sion (few MK) is only an indirect means of inferring outflows 
in X-rays and a signature that is easily confused with contribu- 
tions from accretion. Direct detection of X-ray emission from 
pre-MS jets by means of a displacement with respect to the cen- 
tral coronal sour ce has been achiev ed only in a handful of cases 
(see summary in lBonito et al.ll2007l) . 

In this article we examine the X-ray properties of Z CMa dur- 
ing its recent optical outburst, we compare them to its quiescent 
state, and we present the X-ray detection of its jet. 

ZCMa is a 0.1" pre-MS binary. While the south-east (SE) 
FUor object dominates the light at optical wavelengths, the 
north-west (NW) co mponent is a powerful infrared (IR) source 
dKoreskoet al.|[l991l) . The FUor star has SM© and is lo- 
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Fig. 1. Lightcurve of ZCMa; vertical lines indicate the date of 
Chandra observations. 



Table 1. Observing log for Chandra exposures of ZCMa. 



ObsID 


SeqNo 


Instrument 


Date [UT] 


Expo [ksec] 


3751 


200221 


ACIS-S2 


2003-12-07 22:00 


38 


10845 


200569 


ACIS-S3 


2008-12-28 22:48 


40 



cated just below the birthline in the HR diagram dHartmann et al.l 
11989 ). Assuming that the IR sour ce is coeval, it can be mod - 
eled as a BO III star with 16 Mq dvan den Ancker et al.1 [20041) . 
The optically dominating component is, therefore, the secondary 
in the binary system. This star is likely responsible both for 
the FUor phenomena and for the jet and molecular outflow 
observed at radio and optical wavele ngths (IP oetzel et al.lll989t 
lEvans et al.ll994tlVelazquez & Rodrig uez 200 1]). The long-term 
lightcurve of Z CMa exhibits features of both FUor and EXor- 
like events: Optical outbu rsts of ^ Imag are super posed on a 
two decade long decay (Ivan den Ancker et al ] 12004 . An alter- 
native explanation for the irregular light variations of Z CMa is 
scattered light from the embe dded primary that penetrate s an en- 
velope of variable thickness dHartmann & Kenvonl l 1996V 

In Jan 2008 ZCMa started a large optical outburst (see 
Fig. [U. Visual brightness estimate^ are available for the ini- 
tial gradual increase by ~ 1.5 mag. Further visual observations 
and additional BVRcIc CCD photometric observation^ have 
resumed after a ^ 6 month-long gap. Z CMa reached its (tempo- 
rary) maximum of ~ 8.3 mag at the beginning of Jan 2009. 



2. Data analysis and results 

A Chandra o bservation of ZCMa carried out in Dec 2003 was 
presented by IStelzer et al.l d2006l) in the framework of a sur- 
vey for X-ray emission from Herbig stars. During that observa- 
tion the source was not in optical outburst. We obtained another 
40 ksec of Chandra Director's Discretionary Time to search for 
changes of the X-ray characteristics during the recent optical 
outburst. Table [T] gives the observing log. For consistency we 
re-analyze here the 2003 observation analogous to the new data 
set. 



' We obtained the visual photometry from the American Association 
of Variable Star Observers (AAVSO) at http://www.aavso.org 

^ The CCD photometry was obtained by 
Czech observers L. Brat and L. Smelcer; see 
http://var2.astr o.cz/meduza/light-curves-ccd.php?star=Z^ 20CMa&shv=i 
ofVSES CAS. 
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Fig. 2. Chandra/ ACIS images for the 0.2 — 8keV band in a 
25" X 25" region around ZCMa (0.5" pixels) for the 2003 
(top) and the 2008 (bottom) observation; zoom on Z CMa in the 
0.2 — 1 ke V band for both epochs. All detected X-ray sources are 
marked with black circles, and the binary components are shown 
as white x-shaped symbols. The position of the X-ray source de- 
tected in the 2008 image to the SW of the central source labeled 
'J' is overlaid in the 2003 image as red circle. In 2008, the soft 
band image shows continuous emission along the jet axis rather 
than distinct X-ray sources. 



We used the CIAO packag^E version 4.0, and we started with 
the level 1 events file provided by the Chandra X-ray Center. 
In the process of converting the level 1 events file to a level 2 
events file for each of the observations we performed the follow- 
ing steps: We removed the pixel randomization in order to opti- 
mize the spatial resolution. We filtered the events file for event 
grades (retaining the standard grades 0, 2, 3, 4, and 6), and ap- 
plied the standard good time interval file. 

For our science goal of detecting the X-ray source(s) associ- 
ated with ZCMa, source detection was restricted to a 100 x 100 
pixels wide image (1 pixel = 0.5") and a congruent, monochro- 
matic exposure map for 1.5 keV centered on the optical position 
of ZCMa. Source detection was carr ied out with the WAVDE- 
TECT algorithm (.Freeman et al.ll2002l) . We used wavelet scales 
between 1 and 8 in steps of We tested a range of detection 
significance thresholds, and found that crth = 10~^ avoided spu- 
rious detections and at the same time separated close emission 
components. 

Mci CIAO is made available by the Chandra X-ray Center and can be 
downloaded from http://cxc.harvard.edu/ciao/download/) 
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Table 2. X-ray parameters of the Z CMa binary and its jet. 



Date 


Opt/IR 


Off ax 


Counts 


Expo* 


Pks 






['] 


in 0.3 - 8keV 


[sec] 




Dec 2003 


ZCMa 


1.95 


47.6 ± 6.9 


17857. 


0.54 


Dec 2003 


jet 


1.91 


< 9.4 


17857. 




Dec 2008 


ZCMa 


0.30 


107.2 ± 10.4 


39594. 


0.50 


Dec 2008 


jet 


0.34 


20.7 ±4.6 


39610. 


0.32 



* The effective exposure time in the Dec 2003 observation is 
substantially reduced with respect to the duration of the observation 
because the star is near a chip border. 



The two X-ray images are shown in Fig. |2]with optical po- 
sitions of the binary and detected X-ray sources overlaid. The 
brightest X-ray source coincides in both observations with the 
unresolved binary star. Two new X-ray sources, not seen in the 
2003 data, are detected in the 2008 image. Especially remark- 
able is the faint source in the S W elongation of Z CMa (named 
'J' henceforth). The position angle of this source is (225 ± 5)°, 
roughly in agreement with the orientation of the bl ue-shifted part 
of the optical jet detected by iPoetzel et alj (Il989h . In the broad 
band image (0.2 — 8keV) this source is separated by 2.4" from 
Z CMa. If only photons below 1 ke V are considered the 2008 
image shows elongated emission along the same axis suggesting 
a chain of weak X-ray sources extending to the S W of Z CMa. 
There is no significant soft excess in the SW direction during 
2003 (see zoom in Fig. |2]). The remaining X-ray sources (la- 
bels 'N95-B' and 'N95-C' in Fig. ^ are tentatively identified 
with p re-MS candidates mentioned by Nakaii ma & Golimowskil 
( Il995h . In the following we concentrate on the analysis of the 
X-ray emission associated with ZCMa and its jet. 

We calculated the source count rates in the following way: 
First, the point-spread-function (PSF) was computed for each 
X-ray position. A circular source photon extraction region was 
defined as the area that contains 95 % of the PSF. Only for source 
'J' the extraction region was restricted to 90 % of the PSF to 
avoid contamination from the wings of the central source. The 
background was extracted individually from a squared region 
centered on the source extraction area and several times larger 
than the latter one. Circular areas centered on the positions of the 
X-ray sources were excluded from the background area. The S/N 
was computed from the counts summed in the source and back- 
ground areas, respectively, after applying the appropriate area 
scaling factor to the background counts. In practice, the back- 
ground is very low (a fraction of a count in the source extraction 
area). Finally, count rates were obtained using the exposure time 
at the source position extracted from the exposure map. We have 
estimated a 95 % confidence upper limit for the count rate at the 
position of 'J' in the Dec 200 3 observation (mar ked red in Fig.|2] 
left) using the algorithm of iKraft et all (1199 lb . In Table |2] we 
summarize the relevant X-ray parameters of the Z CMa binary 
and the source 'J' identified with the jet. 

An individual response matrix and auxiliary response were 
extracted for the position of each source using standard CIAO 
tools. The spectrum of the brighter source (ZCMa) was binned 
to a minimum of 10 counts per bin and that of the fainter 
one ('J') to 5 counts per bin. As mentioned above, the back- 
ground of ACIS is negligibly low. We fitted the spectra in the 
XSPEC 12.4.0 environment with a one- or two-temperature ther- 
mal model subject to photo-absorption (WABS • APEC and WABS 
• [APEC + APEC], respectively). 



For the brighter source (ZCMa) the 1-T fit of the Dec 
2008 observation displays substantial residuals slightly below 
IkeV, and we resort to the 2-T model. The best fit A'h is 
Tig • 10^^ cm^^, compatible with the range of va lues pub- 
lished for the optical extinction Ay = 2.4...4.6mag (lE lia et aLl 
l2004tlAcke & van den Anckeij|2004 . Assuming a particle den- 
sity of 0.3 cm~^ the galactic absorption at 1 kpc amounts to 
^ 10^^ cm~^. Therefore, ZCMa may have some additional 
absorption related to the star forming environment. On the 
other hand, our measurement does not rule out negligible cir- 
cumstellar absorption. In our best-fit model the soft compo- 
nent (kTi — OA'^Q^^eV) dominates over the hard compo- 
nent (kT2 ~ 7.5 keV; unconstrained) with an emission measure 
logEMi [cm-3] = 53.9;j| vs. log£;Af2 [cm-^] ^ 53.ltoi- 
Due to the large confidence intervals of the spectral parameters, 
we determine only a lower limit on the X-ray flux, adopting a 
minimum A^h of ^ 10^^ cni^^ coming up for the interstellar 
absorption but neglecting any possible contribution from the en- 
vironment of the star. We find /x > 1.6 • lO"^** erg/cm^/s for 
the 0.5 — 8.0 keV band, corresponding to log [erg/s] > 30.3 
at a distance of 1050 pc. 

For the Dec 2003 observation due to low photon statistics 
we can not formally exclude the 1-T model. Iso-thermal models 
are known to be an oversimplification and the number of ther- 
mal components needed to fit stellar X-ray spectra generally in- 
creases with photon statistics. We do not present a detailed spec- 
tral analysis of this data set because the quality is poor. However, 
we can test the compatibility of the 2003 spectrum with the 2- 
T model parameters derived from the 2008 observation. Indeed, 
we find from a direct comparison of the 2008 best-fit model to 
the spectrum observed in 2003, without fitting it to the data, a 
Xrod ~ 1-0 (5 d.o.f.). The good agreement is also evident from 
Fig. [3] where we plot both observed spectra together with the 
best-fit model (dashed lines) of the 2008 data. Differences in the 
appearance of the two spectra in Fig. [3] are due to the fact that 
the data and the model are shown folded with the instrument re- 
sponse matrix. We recall here that Z CMa is located on two dif- 
ferent CCD chips and at different off-axis angle in the 2003 and 
2008 data sets (see Table [T]|, and this implies different spectral 
response and effective area. 

The fainter source ('J') has only 20 counts in Dec 2008. The 
1-T model gives the same absorption as found for ZCMa al- 
beit with even larger uncertainties (A^h = 7ly • 10^^ cm~^). 
The large error bar of the column density makes it impossible 
to determine the emission measure and flux of this source to 
even an order of magnitude precision. The temperature of source 
'J' is also not well constrained (0.2^°^ keV) but probably rela- 
tively soft. We find a median photon energy of 0.9 keV for 'J' 
versus 1.2 keV for ZCMa. Analogously to the case of ZCMa, 
we estimate a lower limit for the X-ray flux of 'J' adopting A^h 
of ~ 10^^ cm^^ coiTesponding to the expected interstellar ab- 
sorption. The derived value of 2 • 10^^^ crg/cm^/s translates 
to logix [erg/s] > 29.4 if we assume that the distance of this 
source is the same as for the one associated with Z CMa. 

Lightcurves were extracted and searched for variability with 
a maximum likelihood method that divides t he sequence of pho- 
tons in intervals of constant signal (see iStelzer et al.l l2007l) 
and, independently, with the Kolmogorov-Smirnov (KS) test. 
According to this analysis both sources are not variable within 
the individual exposures. The significance level of the KS statis- 
tic for each source is given in the last column of Table |2] where 
the large values of Pks indicate that the data is not significantly 
different from the assumption of a constant source. 
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Fig. 3. Folded Chandra/ ACIS spectra of Z CMa. Overplotted on 
both spectra is the best fit model obtained for the 2008 observa- 
tion (dashed lines). 



3. Discussion 

3. 1. The central X-ray source 

Prior to our observations, only two pre-MS stars have been mon- 
itored in X-rays during an optical EXor outburst. V1647 0ri 
has shown a factor 100 increase of the X-ray luminosity and 
a steady decay in the subsequent two years correlated with the 
near-IR lightcurve. With respect to the pre-outburst state the 
spectrum was hard throughout the decaying phase (3.6 keV) 
and inco mpatible with X-ray emission from accretion shocks 
dKastner et al. 2006) . The other E Xor object observed i n X-rays 
during an outburst is V 1 1 1 8 Ori. lAudard et all (l2005h showed 
that its X-ray luminosity varied only by a factor of two with re- 
spect to the pre-outburst value but the spectrum softened durin g 
the E Xor event. A possible explanation given by Audard et al.l 
(l2005i) was a changing structure of the magnetosphere that 
moved inward as a result of strong accretion and disrupted the 
higher, hotter coronal loops. After the outburst V1118 0ri was 
seen to have faded by a factor of four in X-ray lu minosity but has 
not si gnificantly changed its X-ray temperature dLorenzetti et al.l 
I2006h . 

For Z CMa the X-ray brightness and temperature are remark- 
ably constant when the 2003 quiesence and the 2008 outburst 
data are compared. Overall, the temperature structure of the X- 
ray spectrum of Z CMa is typical of corona l sources observed 
at similar statistics (e.g. lSchmitt et al.lll990h . Clearly, the spec- 
trum is much softer than that of VI 647 Ori during its EXor 
event. The low statistics make it impossible to single out even- 
tual contributions fro m accretion and/or j et shocks. Similarly, 
low statistics impeded lSkinner et al. I (120091) to constrain the tem- 
perature of the soft component in the X-ray spectrum of the 
FUor star V1735Cyg that they observed outside outburst with 
XMM-Newton. Their confidence level on the temperature of the 
soft component of V1735Cyg is similar to ours (0.1 — IkeV) 
leaving open the question on whether this emission is accretion- 
related. We have estimated a lower limit to the X-ray luminosity 
of Z CMa that is compatible with the X-ray luminosities reported 
by lSkinner et all (12009) for V 1735 Cyg and FUOri. As those au- 
thors pointed out, these luminosities are at the high end of the 
distribution of T Tauri stars possibly implying that FUor objects 
have higher than average T Tauri mass. For Z CMa this observa- 
tion is consistent with the high mass of 3 Mq from the literature. 

The binary nature of Z CMa leaves room for further specu- 
lations on the origin of the X-ray emission: (i) The X-rays could 



be from the less massive star, which is in this case non-variable 
both in the optical and X-rays, and the optical variations due to 
changes in scattered light from the primary. However, the recent 
optical photometry has a clear outburst signature and we discard 
the scattered-light scenario, (ii) The optical variations may come 
from the FUor star and the X-rays from the embedded primary. 
This hypothesis is also unlikely as much higher X-ray absorption 
would be expected than is observed. We conclude that probably 
the same star (the less massive optical component) is respon- 
sible for both the optical outburst and the X-ray emission but 
that both phenomena are unrelated. Although the existence of 
a corona can not be taken for granted for an intermediate-mass 
star, our recent detection of the magnetic field of Z CMa supports 
this interpretation (Hubrig et al., 2009, A&A submitted). 

3.2. The displaced X-ray source 

A new X-ray source is detected in the 2008 image in the SW 
elongation of the source associated with the Z CMa binary. Its 
position angle agrees with that of the optical jet (Poetz el et al.l 
[1989) and the radio jet (Velazquez & Rodriguez 2001) suggest- 
ing an interpretation as emission from a n internal jet shock. 
According to the strong shock condition (IZel'dovich & Raizerl 
19661), our confidence interval of measured X-ray tempera- 
tures, 0.1...1keV, is compatible with shock velocities between 
^ 300...900km/s. At a distance of ^ 2", corresponding approx- 
imately to the location of the new X-ray source, in op tical forbid- 
den lines the jet was seen to move at < 600km/s (P oetzel et al.l 
1989). Considering the unknown inclination this is a lower limit 
for the true velocity of the jet propagation. Our estimate for the 
shock speed is, therefore, not incompatible with the expecta- 
tion that Wshock < Vjct- ZCMa has a higher mass than most 
other pre-MS stars with X-ray detected jets (see summary by 
Bonito et al. 2007). The large distance of the displaced X-ray 
emission from the central star is remarkably similar to that of 
the only o ther high-mass pre -MS object with known X-ray jet, 
HH 80/8 1 dPravdo et al.l2()04l) . For the X-ray luminosity we have 
estimated a conservative lower limit of 2.6 • 10^^ erg/s, that could 
easily be underestimating the intrinsic luminosity by one order 
of magnitude or more if the extinction was much higher than 
the assumed interstellar value. Our result is compa tible with X- 
ray luminosities derived for other pre-MS jets (see lBonito et al.l 

2007) . However, all X-ray observations of this object class are 
hampered by similarly low statistics and, consequently, poorly 
constrained parameters. 

Due to the relatively low X-ray temperature that can be 
achieved in pre-MS jets an image limited to soft photons should 
better trace the outflow. Indeed, we find continuous soft emis- 
sion between ZCMa and source 'J' in the 2008 image. Possibly 
several knots along the jet axis are producing X-rays at a time. A 
similar finding was reported for DG Tau where the X-ray jet was 
detected both at close (50 AU) and larger (1140AU) separation 
from the central star dGiidel et al. 1 120081; [Schneider & Schmit3 

2008) . 

In 2003 there was no X-ray source at the position of ' J' . The 
upper limit to the count rate at this epoch does not exclude X- 
ray emission at similar levels as observed in the 2008 image. 
However, examination of the photon centroids reveals that all 
counts are in the upper left of the photon extraction area (see 
Fig. |2|l indicating weak emission in the same SW direction as 
'J' but at a distance of only ^ 1.2" from the central source. If 
this emission is interpreted as the head of a moving shock front 
that has expanded (and brightened) within the last 5 years to 
the position of 'J', its projected velocity would be ~ 9GGkm/s. 
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This value is in rough agreement with our measurement from 
the X-ray temperature. However, the analysis of the soft-band 
image from 2008, where continuously extended soft emission 
from ZCMa out to ^ 2.4" is detected, suggests that 'J' repre- 
sents the temporary radiative loss of a localized blob of material 
impacted by an intermittent jet. The weak emission closer to the 
star may represent similar but independent events. 

We add that alternative explanations are perceivable for the 
X-ray emission displaced from Z CMa: Scattered ste llar light re- 
flected from the outflow cavities was discussed by iBallv et aU 
(I2003h for the case of HH 154 but seems unlikely to hold for 
Z CMa because of the large displacement o f > 2000 AU with re- 
spect to the central star Again for HH 154.l Murphv et all (l2008l) 
suggested magnetic reconnection in the space between two in- 
teracting jets from the two binary components. There is no ob- 
served evidence for a second jet in Z CMa. High-quality high- 
resolution imaging of the Z CMa outflow should help to clarify 
this question. Monitoring of the X-ray morphology during the 
next decade should also help to constrain the nature of the emis- 
sion. For an assumed jet diameter of 100 AU we can estimate the 
electron density from the observed X-ray emission measure and 
derive a radiative cooling time of ^ 20 yrs. 
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